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AWARDS ABSTRACT 

The invention relates to a method for providing a polarization 
filter for processing synthetic aperture radar (SAR) image data, and more 
particularly to a method for maximization of signal-to-noise ratio and 
discrimination between areas of different target types by utilizing the 
vector nature of the return waves in order that information about the 
surface contained in the polarization properties of the scattered return 
waves may be recovered, 

A systematic theoretical analysis and development of a polarization 
filter is described for two applications: (1) maximization of signal-to- 

noise ratio and (2) discrimination between two target types. The method 
employed finds the optimum receive antenna polarization analytically, u 
relies on a numerical solution to find the optimum transmit antenna polari- 
zation. The analysis uses the Stokes matrix representation, and therefore 
can be employed to analyze the partially polarized scattered field from 
extended targets and can also be applied to the bistatic case. Examples o 
the technique are presented for the monostatic NASA CV990 polarimetric 
L-band radar. FIG. la illustrates the total SAR power image (also known as 
the Span of the Stokes matrix). FIG. 1b illustrates the outlines of dif- 
ferent areas of target types and the optimum filter image of the SAR power 
image of FIG. la. FIG. 1c illustrates the enhancement factor image which 
is the ratio of the two other images. FIG. Id shows the outlines of dif- 
ferent areas of target types and gives the position of two 3tu ^ied at. 
No 1 is a forested area site and No. 2 is an urban area site. FIGs. 2 ana 
3 illustrate an optimization signature maximization of signal-to-noise 
ratio for an urban area and a forested area, respectively, and F . 
illustrates an optimization signature maximization of contrast ratio be- 
tween urban and forested areas. FIG. 5 is a graph of enha . n 7f® nt 1 Q f Q a Q Ct °£ 
variation with range line. The line 0 is the near range and line 1000 is 
the far range of the SAR image. FIG. 6 is a graph of enhancement factor 
variation with size of training area. The training area contains 
Dixels FIG. 7a illustrates the total power image (also known as Span t 
the Stokes matrix) of the Blackwater River area. FIG. 7b is the optimum 
filter image of FIG. 7a, and FIG. 7c is the enhancement factor image of 
FIG. 7a, which is the ratio of the two other images. FIG. 8 illustrates an 
optimization signature maximization of signal-to-noise ratio for an urban 
area with symmetrized noise. 

The novelty of the invention resides in a particular method of 
polarization filtering based on Stokes matrix representation to enhance 
certain characteristics of SAR images. 
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METHOD FOR PROVIDING A POLARIZATION FILTER FOR 
PROCESSING SYNTHETIC APERTURE RADAR IMAGE DATA 


Origin of the Invention 

The invention described herein was made in the perform- 
ance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
Contractor has elected not to retain title. 


Technical Field 

The invention relates to a method for providing a 
polarization filter for processing synthetic aperture radar 
(SAR) image data, and more particularly for a filter for 
maximization of signal-to-noise ratio and a filter for dis- 
crimination between areas of different target types by utiliz- 
ing the vector nature of the return waves in order that 
information about the surface contained in the polarization 
properties of the scattered return waves may be recovered. 

Background Art 

Conventional imaging radars operate with a single 
fixed-polarization antenna for the radio frequency signal 
transmission and reception. As a consequence, for every 
resolution element (pixel) in the image, the scattered wave (a 
vector quantity) is measured as a scalar quantity and any 
additional information about the surface contained in the 
polarization properties of the reflected signal is lost. To 
insure that all the information in the reflected signal is 
retained, the reflected signal must be measured as a vector, 
i.e., both the amplitude and the phase should be measured. 
The greater information derived from the polarized measure- 
ments provides a more complete description of the backscatter 
phenomena of the target area. This greater information can, 
for example, be used in target discrimination, target classi- 
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fication or feature enhancement. This invention relates to 
techniques for utilizing polarization information to enhance 
certain characteristics in SAR images. 

Determination of the optimal polarization state to 
5 enhance an image has been extensively studied in the past. The 
scattering matrix co-polarization and cross-polarization nulls 
represent a solution when the returns from a single point 
target need to be minimized. A. B. Kostinski and W. M. 
Boerner , "On the Polarimetric Contrast Optimization," IEEE 
10 Trans. Antennas and Propagation, Vol. AP-35, No. 8 , pp. 
988-991, August 1987, based their analysis on the Graves power 
matrix to determine the optimum transmit and receive polariza- 
tions. This technique is only applicable to maximize the 
contrast between two specified point targets and has been 
15 applied to polarimetric radar images by estimating an equiva- 
lent scattering matrix representation for an extended area in 
an image. B. James, A. B. Kostinski and W. M. Boerner, "Po- 
larimetric Matched Filter for POLSAR Image Interpretation of 
Ocean Surface Scatter," Proc. IGARSS ' 88 , P. 67, Edinburgh, 
20 United Kingdom. 

The accuracy of this technique for extended targets is 
unknown since representation requires either an average Stokes 
matrix [J. J. van Zyl, A. Zebker and C. Elachi, "Imaging radar 
polarization signatures: Theory and observation," Radio Sci- 
25 ence, 22(10, pp. 529-5^3 , July/Augu 3 t 1987; and J. J. van Zyl, 
C. H. Papas and C. Elachi, "On the optimum polarizations of 
incoherently reflected waves," IEEE Trans, on Antennas and 
Propagation," Vol. AP-35, No. 7, July 1987] or an average 
covariance matrix [J. A. Kong, A. A. Swartz, H. A. Yueh, L. M. 
30 Novak and R. T. Shin, "Identification of Terrain Cover Using 
the Optimum Polarimeter Classifier," J. Electromagnetic Waves 
and Applications, Vol. 2, No. 2, pp. 171-1 9 1 !, 1988]. Both 

these representations use the second order statistics of the 
scattering matrix. G. A. Ioannidis and D. E. Hammers, "Opti- 
35 mum Antenna Polarizations for Target Discrimination in Clut- 
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ter," IEEE Trans. Antennas Propagat., Vol. AP-27, No. 3, May 
1979, introduced a method based on Lagrangian multipliers to 
solve for the optimal polarization using the Stokes matrix. 
However, for some cases their solutions violate the constraint 
that the Stokes vector for the receive antenna must be fully 
polarized. Additionally, A. A. Swartz, H. A. Yueh, J. A. Kong, 
L. M. Novak and R. T. Shin, "Optimal Polarizations for achiev- 
ing Maximum Contrast in Radar Images," J. Geophys. Res., Vol. 
93, No. B12, pp. 15252-15260, December 1988, developed a 
parallel method based on the covariance matrix. However their 
analysis is restricted to the backscatter case. 

Statement of the Invention 

An object of this invention is to provide a polariza- 
tion filter for maximization of the signal-to-noise ratio of a 
single target return or maximization contrast between two 
target types, particularly between man made urban areas and 
natural nonurban areas (forest, grass, and ocean). 

In accordance with the present invention, polarization 
filters based on Stokes matrix/Stokes vector representation 
maximizing the signal-to-noise ratio are developed for differ- 
ent noise characteristics and maximizing the contrast between 
different target types by maximization of the power ratio 
between target types. 

Brief Description of the Drawings 

FIG. la illustrates the total SAR power image (also 
known as the Span of the Stokes matrix). FIG. 1b illustrates 
the outlines of different areas of target types and the opti- 
mum filter image of the SAR power image of FIG. la. FIG. 1c 
illustrates the enhancement factor image which is the ratio of 
the two other images. FIG. Id shows the outlines of different 
areas of target types and gives the position of two studied 
sites; No. 1 is a forested area site and No. 2 is an urban 
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FIG. 2 illustrates an optimization signature maximiza- 
tion of signal- to- noise ratio for an urban area. 

FIG. 3 illustrates an optimization signature maximiza- 
tion of signal-to-noise ratio for a forested area. 

FIG. 1 ! illustrates an optimization signature maximiza- 
tion of contrast ratio between urban and forested areas. 

FIG. 5 is a graph of enhancement factor variation with 
range line. The line 0 is the near range and line 1000 is the 
far range of the SAR image . 

FIG. 6 is a graph of enhancement factor variation with 

p 

size of training area. The training area contains n pixels. 

FIG. 7a illustrates the total power image (also known 
a 3 Span of the Stokes matrix) of the Blackwater River area. 
FIG. 7b is the optimum filter image of FIG. 7a, and FIG. 7c is 
the enhancement factor image of FIG. 7a, which is the ratio of 
the two other images. 

FIG. 8 illustrates an optimization signature maximiza- 
tion of signal-to-noise ratio for an urban area with symmet- 
rized noise. 

Detailed Description of the Invention 

A systematic theoretical analysis and development of a 
polarization filter is described for two applications: (1) 
maximization of signal-to-noise ratio and (2) discrimination 
between two target types. The method employed finds the 
optimum receive antenna polarization analytically, but relies 
on a numerical solution to find the optimum transmit antenna 
polarization. The analysis uses the Stokes matrix representa- 
tion, and therefore can be employed to analyze the partially 
polarized scattered field from extended targets and can also 
be applied to the bistatic case. 

Examples of the technique are presented for the mono- 
static NASA CV990 polarimetric L-band radar which alternately 
transmits illuminating waves with horizontal (h) and vertical 
(v) polarization while the receiver employs two channels, one 
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for the horizontally polarized scattered waves (h' ) and the 
other for the vertically polarized scattered waves (v') in 
order that the data received and stored may be processed in a 
bistatic scattering matrix 
5 

s h'h s h'v 
^v'h ^v'v 

where: S h , h is the complex ratio of the electric field of the 
10 horizontally polarized part of the scattered wave ( h T ) and the 
horizontally polarized part of the illuminating wave (h); 

S^, v is the complex ratio of the electric field of the 
horizontally polarized part of the scattered wave ( h T ) and the 
vertically polarized part of the illuminating wave (v); 

15 Sy'h is the complex ratio of the electric field of the 

vertically polarized part of the scattered wave (v*) and the 
horizontally polarized part of the illuminating wave (h); and 
Syt v is the complex ratio of the electric field of the 
vertically polarized part of the scattered wave (v’) and the 
20 vertically polarized part of the illuminating wave (v). 

The terminology and definitions used here and herein- 
after have previously been defined by J. J. van Zyl, H. A. 
Zebker and C. Elachi, (July/August 1987) and are Incorporated 
herein by reference. 

25 Image enhancement filters maximizing the signal-to- 

noise ratio are developed for different noise characteristics 
and different target types (urban and forest). A filter is 
also developed to maximize the power ratio between urban and 
natural (forest or ocean) targets for optimum discrimination 
30 between the types of targets. Results are presented to demon- 
strate that the filter maximizing the contrast between urban 
and ocean areas i3 essentially the same as the one maximizing 
the contrast between urban and forested areas. 
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Matched Filter In Presence of Background Noise 

Let H be a scattering matrix characteristic of the 
received noise. 


N h'h N h'v (i) 

Ny'h Ny’v 


Each element of N is a complex random variable with zero mean 
and a known variance a 2 (the noise power). Assume that the 
elements of N are uncorrelated and that the noi3e of each 
10 radar channel has the same power. The mathematical formula- 
tion is as follows: 

<N ij N kl >= ° 2 lf i=k and J =1 
-0 otherwise 

<N^j>=0 for all i and j 

15 where < > denotes the spatial average. 

The corresponding average Stokes matrix is of the form: 


20 


"n 


o 2 0 0 0 

0 0 0 0 
0 0 0 0 
0 0 0 0 


( 2 ) 


Equation (2) states that the average noise signal is unpolar- 
ized when the noise is uncorrelated. Therefore, the noise 
merely contributes a constant power level to the signal [van 
Zyl , et al., (July 1987)]. 

25 To find the polarization configuration of the receive 

and transmit antenna maximizing the signal-to-noise ratio 
(SNR), let A be the average Stokes matrix characterizing a 
target of interest, and "n be the average Stokes matrix for 

A A 

the noise. Also let S p and S*. be the Stokes vectors charac- 
30 terizing some arbitrary polarization configuration of the 
receive and transmit antennas. 


s to 


s r0 

s*. 

and S r * 
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where 


( 3 ) 


Vi 

3 r2 

3 r3 


and s. 


Vi 

3 t2 

3 t3 


Because both Stokes vectors represent an antenna configurati- 
on, a necessary condition [van Zyl, et al. t (July 1987)] is: 


5 t0 


,2 

5 r0 


3 r * 3 r ’ 


(4) 


10 


15 


20 


25 


It can be shown [van Zyl, et al., (July 1987)] that the sig- 
nal-to-noise ratio SNR is given by: 




A S, 


SNR 


^ rn A 

s r N s t 


(5) 


Without loss of generality, it can be assumed that S r Q-S^.Q=1 
since these terms cancel in the ratio. 

The received power for a given antenna polarization can 
be maximized or minimized by maximizing or minimizing S p q and 
S to . However, in practice S r0 and S tQ are proportional to the 
gain of the antennas, which for a given wavelength is deter- 
mined by their sizes and cannot easily be varied. Consequent- 
ly, the present invention uses signal polarization 
characteristics to extract information from SAR data assuming 
that the gain of the transmit and receive antennas are fixed. 
Therefore, Equation (4) can be rewritten as: 


s p • s p = 1 and s t • s t - 1 (6) 

Assuming uncorrelated noise, the denominator of Equa- 
30 tion (5) is constant and equal to o^. Therefore, maximizing 
the signal-to-noise ratio is equivalent to maximizing the 
numerator. To optimize SNR, the average Stokes matrix A is 
written as follows: 
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5 


10 


15 


A = 


m u 

A 

v Q 


(7) 


where u denotes the transpose of u and u, v are 3-element real 
vectors and Q is a 3*3 real matrix. 

Using Equation (7) for A, an alternate form for the 
numerator of Equation (5) is: 


s£ A S t - m+u-s t +v*s r +s r *Qs t (8) 

This is the expression to be maximized under the constraints 
expressed in Equation (6), namely the assumption that the gain 
of the transmit and receive antennas are fixed. However, the 
solution of constrained maxima for any one of the variables 
may be extremely difficult. It is therefore desirable to use 
the Lagrangian multiplier method to solve this problem. After 
the necessary differentiation, the s^ and s p that maximizes 
the SNR is the solution of the equations: 


20 


G+Qs r - A 1 s t 

(9) 

v+Qs t - X 2 s p 

(10) 


where A 1 and X 2 are the Lagrangian multipliers. 

If the transmit antenna is known (i.e., s t is given), 
the receive antenna polarization which maximizes the SNR is 
determined from Equations (6) and (9) to be as follows: 


r 

I|V * 03 t ll 

where || || denotes the norm. Physically, Equation (11) states 

that the power is maximized when the polarization of the 
receive antenna is matched to the polarization of the scat- 
30 tered wave incident upon the receive antenna. 

The transmit antenna polarization s^. must then satisfy: 

(Q Q — Ai X 2 I)Sj. * v — ^2^* (12) 
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When Equation (11) is satisfied, the numerator of Equation (5) 
can then be written: 

p m ax * m + u-s t + ||v+Qs t || (13) 

5 where P__„ is the maximum obtainable scattered power given the 

ulaX 

transmit antenna polarization. 

It is hard to analytically solve for the transmit 
polarization which would maximize P max * Instead, this optimum 
transmit polarization may be found numerically, as will be 
10 shown below under the heading Numerical Examples. 

Matched Filter as a Discriminator Between Two Target Classes 

One objective of this invention is to find the polari- 
zation of the receive and transmit antennas which maximizes 
15 the ratio of signal power scattered by one type of target to 

that scattered by another type of target. Let F^ represent 

( 2) 

the Stokes matrix characteristic of target No. 1 and F 
represent the Stokes matrix characteristic of target No. 2. 
The problem then Is to find the Stokes vectors S r and S^ 
20 representing the receive and transmit antenna polarization 
realizations which maximize: 


p (1) s£ F (1) S t 

— 1 (14) 

25 P (2) F (2) S t 

But first it is necessary to solve the following problem. 
Assuming the transmit antenna polarization is fixed, only the 
receive antenna is to be operated upon to maximize the power 

A A 

ratio given by Equation (14). S 1 and S 2 are defined as 
30 follows: 


( 15 ) 
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F<«i t 


5 02 


(16) 


S r , the Stokes vector of the receive antenna, can be written 
as: 




(17) 


As before, since S r represents the Stokes matrix of an an- 
tenna, it is a completely polarized wave and hence must sat- 
isfy: 


Sp 0 - Sp*s r - 1 (18) 

It is clear that to maximize the ratio given by Equation (1*0, 
it is necessary to maximize the contrast 


C (s r ) 


s 01 +3 r* 3 1 

S 02 +3 r* 3 2 


(19) 


Since C(s r ) is the ratio of two received powers, it is reason- 
able to insist that: 


C (Sp) * 0. (20) 

To maximize the contrast between the two targets, one has to 
maximize the return from the one target, while at the same 
time minimizing the return from the other target. For the 
moment, assume that there is no antenna polarization which 
could cause the denominator of Equation (19) to become zero 
and then generalize the result. To optimize C(s r ) subject to 
the constraint condition given by Equation (18), Lagrangian 
multipliers are introduced. After performing the necessary 
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differentiations, it is found that the optimum polarizations 
are the solutions to 


5 


( s oi + V 3 i\ 

3l " 32 \ S 02 +s r* 3 ^ 
where y is the Lagrangian multiplier, 
term 


y3 p , (21) 

Introducing a shorthand 


S 0 i + 3 r * s 1 

10 a - , ( 22 ) 

S 02 + 3 r * 3 2 

Equation (21) can be rewritten to read 

s(a) - us r (23) 

15 

where s(a) z s^- aSg. 

From Equations (23) and (18) it follows that 
20 y - ±/s(a)‘s(a), ( 2 H) 


and the optimum polarizations are therefore described by 

3 (a) 

3 r * i 7 7 V (25) 

25 r /s(a)*s(a) 

Notice that s r , as expressed in Equation (25), always repre- 
sents a fully polarized wave consistent with the constraint in 
Equation (18). A previous analysis by Ionnadis and Hammers 

(1979) found the optimum receive antenna Stokes vector: 

30 

S r - - R(F ( 1 ) -pF ( 2 ) )S t - - S.-^Sp (26) 

1 y * u y 1 y c 

A fs 

where p Is a Lagrangian multiplier. In general, S 1 and S 2 
35 represent partially polarized waves If ^ and 


are 
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Stoke 3 matrices of extended targets. There is therefore no 

A 

guarantee that the difference S p would be a fully polarized 
Stokes vector. 

In Equation (25), a is still unknown. That unknown can 
5 be solved by substituting the expression for s p in Equation 
(25) into the definition of a given in Equation (22). This 
produces the optimum values of a that are the roots of the 
quadratic: 

10 ( Sq2 _ S 2 * 32 ) 2 ( Sqi Sq 2 _ s ^ * 32^ ^ ^01 — ^ 1 * ^ 1 ^ ® (27) 

If Sq 2*S2*S2, the optimum values of a are 


15 

a max 


S 01 S 02 -s 1 * s 2 


S 02 “ 3 2* 3 2 



^01 S 02 -s 1 * s 2l 


S 02 - 3 2* s 2 


s or 3 r 3 i 


s 02 " 3 2* 3 2 


(28a) 


and 


20 


25 


30 


min 


S 01 S 02~ 3 1 * s 2 
S 02 “ 3 2* s 2 


1 


S 01 S 02~ s 1 * s 2 1 


S 02 - 3 2* 3 2 


s or 3 r 3 i 


S 02 _ 3 2* 3 2 


(28b) 


From Equations (19) and (20), it is clear that a is also the 
ratio to optimize, and therefore Equations (28a and b) also 
define the optimum values of the contrast ratio, C(s p ). It 
can be proven that the contrast ratios given in Equations (28) 
are real, positive and definite values. 

Substituting Equations (28a) and (28b) in Equation (23) 
and substituting the resulting expression in Equation (25), it 
is found that only the solutions with positive signs lead to 
self-consistent results for the contrast ratio. Therefore, the 
optimum antenna polarizations are: 


3 ^ a max ) 

3 


35 


rmax 


ll 3 ^ a max^ll 


(29a) 


F89/411 


13 


and 


(29b) 


3 ^ a min^ 

° r ” ln ' iNwi' 

Now notice that when Sq 2 =s 2 *s 2 , the scattered wave in 
the denominator of Equation (19) is completely polarized. The 
maximum contrast in this case is obviously obtained when an 
antenna polarization which nulls the received power from the 
completely polarized wave is used to receive the scattered 
waves, leading to an infinite contrast ratio. The receive 
antenna polarization is thus chosen to be orthogonal to the 
completely polarized scattered wave: 


s 2 

3 ™ax " ” /s 2 .s 2 


(30) 


To find the minimum contrast ratio for this case, return to 
Equation (23) and find that the quadratic reduces to the 
first order polynomial with 


3 or 3 r 3 i 

ct , « — (3D 

min 2(S 01 S 02 -s r s 2 ) 

Notice that « m j n as given by Equation (31) is the ratio of two 
positive quantities and therefore satisfies Equation (18). 

Note that if the scattered wave in the numerator of 
Equation (18) is completely polarized, the minimum contrast 
ratio is zero which is achieved by using the antenna polariza- 
tion which nulls the received power from this completely 
polarized wave. When both waves are completely polarized, the 
minimum contrast (which is zero) is achieved by using the 
antenna polarization which nulls the received power from the 
completely polarized wave in the numerator of the contrast 
ratio. The maximum contrast (which is finite) is achieved by 
using the antenna polarization which nulls the received power 
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from the completely polarized wave in the denominator of the 
contrast ratio. 

Finally, it should be noted that if the F^' matrix 
appearing in the denominator of Equation (14) is a noise 
5 matrix as defined in Equation (2), then from Equations (16) 
and (21) S 2 -O and 



(32) 


10 where s^v+Qs,,. This result agrees with Equation (2). 

The equations derived in this section describe the 
optimum receive antenna polarization for maximum contrast 
ratio given the transmit antenna polarization. Furthermore, 
Equation (28a) expresses the maximum contrast ratio obtainable 
15 given the transmit antenna polarization state. The next 
logical step is to vary the transmit antenna polarization to 
maximize the contrast over all possible configurations. 
Attempts to find an analytical solution for the optimum trans- 
mit polarization lead to extremely complicated nonlinear 
20 equations. Instead of solving these, this invention presents 
a simpler hybrid method in which the optimum transmit polari- 
zations are found numerically by calculating the maximum 
and/or minimum contrast ratio. The corresponding receive 
polarizations are calculated using Equations (28) and (29) for 
25 all possible transmit polarizations. This hybrid method is 
illustrated with several examples in the following section. 


Numerical Examples 

The procedures developed in the previous sections were 
30 tested with data acquired with the NASA CV990 SAR operating at 
L band. This radar is monostatic and emits alternately hori- 
zontally and vertically polarized waves and receives both 
returns with two co-located antennas, one horizontally polar- 
ized and the other vertically polarized. The system is there- 



F89/411 


15 


fore able to measure a full scattering matrix for each of the 
12x4 meter resolution cells. 

The first test was to minimize the impact of background 
noise using the matched filter approach in the presence of 
5 background noise described above. A second test was to demon- 
strate the use of a polarimetric matched filter to discrimi- 

nate between man-made targets and natural targets. A SAR 
image shown in FIG. 1 of San Francisco acquired in 1985 was 
used in both examples. It was chosen because of the diversity 
10 of the target types as well as the availability of the ground 

truth. FIG. la is the total power image, i.e., the (1,1) 

element of the Stokes matrix. FIG. 1b is the matched filter 
output simulating 45° linearly cross- polarized radar. FIG. 1c 
is the enhancement factor image given by Equation (34) below. 
15 FIG. Id is an outline of different major target areas, namely 
urban, forested and ocean. The data presented in the follow- 
ing examples have been relatively calibrated as that term is 
defined by Zebker, et al., (1986), but not absolutely cali- 
brated. 

20 

Example A : Signal-to-noise ratio maximization (vegetation, 

man-made targets and open water) 

First the matched filter approach is applied to maxi- 
mize the signal to noise ratio. For a given target, the 
25 optimum receive polarization can be analytically computed for 
each transmit polarization. If implemented with sufficiently 
small increments, a numerical search over all polarizations of 
the transmit antenna will yield a global maximum of the sig- 
nal-to-noise ratio as well as the corresponding receive an- 
30 tenna polarization configuration. This search process can be 
visualized through an "optimization signature" which is here 
introduced and defined as a 30 plot where the two horizontal 
axes describe the polarization of the transmit wave, one axis 
being the ellipticity angle, the other the orientation angle. 
35 The maximum signal-to-noise ratio, calculated using Equation 
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(12) at each transmit polarization is displayed along the 
vertical axis. The highest point of the surface will there- 
fore be the global maximum of the signal-to-noise ratio while 
its horizontal coordinates describe the polarization of the 
optimal transmit antenna. The receive antenna polarization is 
computed with Equation (9). 

An optimization signature is displayed for two differ- 
ent types of targets. For the man-made target area No. 2, 
shown in FIG. Id, the Stokes matrix is given by: 


291 .94 

17.35 

29.42 

14.66 

17.35 

229.70 

83.09 

- 12.73 

29.42 

83-09 

- 42.83 

- 22.89 

14.66 

- 12.73 

- 22.89 

105.07 


The corresponding optimization signature indicates that the 
maximum signal-to-noise ratio occurs for a transmit polariza- 
tion of about 10° orientation angle and 0° ellipticity angle. 
The exact numbers are ip 1 0.0°; x r “0.0°;\|) r -9. 8° ; x r ** _ 5° . 

For the natural target area, the Stokes matrix is: 


63.63 

- 6.83 

4.87 

1.29 

- 6 i 83 

32.17 

1.64 

- 3.82 

4.87 

1.64 

1 .92 

3-33 

1.29 

- 3.82 

3.33 

12.29 


25 The corresponding optimization signature shown in FIG. 3 
indicates that the maximum signal-to-noise ratio occurs for a 
transmit polarization of about 85° orientation angle and 5° 
ellipticity angle. The exact numbers in this case are 
ip t =85 . 0° ; x r '5 . 0° ; i|> r =84 . 7° ; x r — 5 • 6° . 

30 

Example B : Discrimination between man-made and natural areas. 

The matched filter approach described above as a dis- 
criminator between two classes is used a3 follows: Two typi- 

cal target areas are chosen as training sites, the first a 
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forested area No. 1 and the second an urban area No. 2 indi- 
cated in FIG. Id by small rectangles. The two areas are 
similar in both size and radar viewing angle. The correspond- 
ing polarization matched filter is computed using the two 
experimentally determined Stokes matrices. An optimization 
signature, as defined in the previous paragraph, is determined 
and presented in FIG. * 1 . For this test, the vertical axis is 
the normalized contrast ratio defined by Equation (13) where 
the maximum contrast ratio is normalized to 1 . The peak of 
the surface occurs for a transmit polarization with 7.5° 
ellipticity angle and H7.5 0 orientation angle. This procedure 
is repeated for twenty target pairs. The polarization of the 
optimal transmit and receive antennas is shown as a scatter 
plot in FIG. 5. The resulting orientation and ellipticity 
angles cluster well around i|)p =1 15 0 and Xr = 0°» tj = )35° and 
X T =0°. This corresponds to a linearly cross- polarized radar 
oriented at 45° from horizontal. As a measure of the perform- 
ance of the filtering method, we define an enhancement factor 
as follows: 


,urb 


Si F 


urb 


P 


for 


t r 1 1 


(35) 


pfor qT pfor n pur b 

P S r F s t P 11 


The enhancement factor is therefore a measure of how radar 
polarization provides discrimination between targets as a 
standalone tool (i.e., without utilizing the target character- 
istics ) . 

FIG. 5 shows that for this example the enhancement 
factor is about 3.0 dB and is range independent. The size of 
the training area does not seem to have an influence on the 
enhancement factor as shown in FIG. 6. This provides a quali- 
tative representation of polarization image enhancement. 

This procedure was also applied to maximize the con- 
trast between urban and ocean areas. The resulting optimum 
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filter is found to be essentially the same as the one optimiz- 
ing the contrast between urban and forested areas. FIG. 7 
displays the effect of the filter on an ocean scene. FIG. 7a 
is the total power image, FIG. 7b is the filtered image simu- 
lating a 45° linearly cross-polarized radar and FIG. 7c is the 
enhancement factor. This filter seems to be a good 
urban/natural discriminator independent of the type of natural 
targets . 

If the uncorrelated noise described with reference to 
Equations (1) through (12) is symmetrized before averaging, 
the corresponding Stokes matrix will have a different form. 
Let N be the symmetrized scattering matrix. We have: 


15 





N »h> 


(36) 


Therefore, the corresponding Stokes matrix, Mjj is written: 


20 



3a 2 0 0 0 

0 a 2 0 0 

0 0 a 2 0 

0 0 0 a 2 


(37) 


Maximizing the signal-to-noise ratio in this case results in 
25 implementing the procedure developed with reference to Equa- 
tions (14) through (32) with being the signal Stokes 

matrix of (e.g., a man-made target) and F v ’ being Mjj. A test 
was conducted on the man-made target Stokes matrix F^ u ^ de- 
scribed with reference to Equation (33) in the previous para- 
30 graph. The optimization signature FIG. 8 is to be compared 
with the optimization signature of FIG. 2. The effect of 
symmetrized noise is clearly noticeable in FIG. 8. The opti- 
mum polarization configuration for the receive antenna and the 
transmit antenna is now i^-l 47.5° , 9° , 

X r -3.8°. 


35 
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Conclusion 

The most Important result of this Invention is the 
development of an optimum polarization filter either for (1) 
maximization of a single target return or (2) contrast en- 
5 hancement between two target types. Because the filter is 
based on the Stokes matrix representation, it is not re- 
stricted to point targets and can be applied to extended 
targets. The filter is also equally valid for bistatic and 
monostatic radar data. 

10 The invention is first embodied in a filter which 

maximizes contrast between urban and natural targets. The 
same filter was found to be optimum when the natural target is 
water, forest, or grass field. Results show that this optimum 
filter improves the contrast ratio by about 3 dB, and that 
15 this enhancement is essentially independent of the incidence 

angle. 

The optimization signature discussed is a useful tool 
as it indicates both the variation and maxima of the contrast 
ratio as a function of the receive antenna polarization. The 
20 polarization resulting in the maximum contrast ratio defines 
the optimum filter. The robustness of this filter is related 
to the variation of the optimization signature surface in the 
neighborhood of the maximum contrast polarization. For exam- 
ple, a rapidly varying surface implies that the filter per- 
25 formance degrades quickly when small errors occur in 

polarization realization. Such small errors can result from 
system imperfections such as crosstalk. 
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METHOD FOR PROVIDING A POLARIZATION FILTER FOR 
PROCESSING SYNTHETIC APERTURE RADAR IMAGE DATA 


ABSTRACT OF THE DISCLOSURE 

Apewrt '-coXcr 

A polarization filter can maximize the signal-to-noise 
ratio of a polarimetric.fSAR' and help discriminate between 
targets or enhance image features, e.g., enhance contrast 
between different types of target. The method disclosed is 
based on the Stokes matrix/ Stokes vector representation, so 
the targets of interest can be extended targets, and the 
method can also be applied to the case of bistatic polari- 


metric radars. 
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